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Abstract 

The time-dependent variational principle proposed by Balian and Veneroni is 
used to provide the best approximation to the generating functional for multi-time 
Green's functions of a set of (bosonic) observables Q M . By suitably restricting the 
trial spaces, the computation of the two-time Green's function, obtained by a second 
order expansion in the sources, is considerably simplified. This leads to a tractable 
formalism suited to quantum fields out of equilibrium. We propose an illustration 
on the finite temperature $ 4 -theory in curved space and coupled to gravity. 
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1 Introduction 



In the past few years, a time-dependent variational principle was proposed M in 
order to deal with the question of how to predict correctly, within an approximation 
scheme, the expectation value of a measured observable, when the system is in a mixing of 
states and is therefore described by a density operator. The variational principle has been 
applied to a variety of quantum processes including heavy ion react ions 0], quantum fields 
out of equilibrium!!, attempts to go beyond the gaussian approximation for fermion|4[] 
and boson systems as well as studies of tunneling through a gaussian barrier [[J . 

More recently, a more elaborate version of the variational principle was developed 
in order to evaluate the generating functional for multi-time correlation functions in equi- 
librium and non-equilibrium statistical mechanics. 

The purpose of the present paper is to extend the formalism of || to quantum 
field theory and to make it more transparent in this area. Among the net advantages of 
the approach proposed here and in [H], is the ability to handle in a consistent and non- 
perturbative way out-of-equilibrium phenomena such as the quantum non-equilibrium 
evolution of the inflaton field(s) in the inflationary scenario^] or that of the quark-gluon 
plasma g. ||. 

We will concentrate in this paper on boson fields because of their direct connections 
with the previously mentioned problems. Extension to fermion fields is quite trivial. 

The paper is organized as follows: in section 2, we recall some results concerning the 
exact determination of multi-time correlation functions of a set of observables. In section 
3, we present the variational action-like which leads to the desired generating functional 
as its stationary value. A peculiarity of the variational principle is the appearance of 
two variational objects akin to a density operator and to an observable. This is however 
analogous to the bra and ket in ordinary quantum mechanics or to the position and its 
conjugate momentum in classical dynamics. In section 4, we select for the variational ob- 
jects gaussian operators and derive the corresponding evolution equations. Being coupled, 
these equations are quite complicated because of the mixed boundary conditions inherent 
in the variational principle. Fortunately, the computation of the one and two-time corre- 
lation functions does not require the solution of the whole set of equations. Indeed, we 
show in section 5 that it is sufficient to expand these equations up to second order in a 
set of sources, and this transforms the problem into two initial value problems. 

In order to illustrate the formalism, we choose in section 6 to focus on a 0(1) 
quantum field theory in Robertson- Walker metric. Furthermore, we restrict the set of 
relevant observables (that is, those we wish to calculate the correlation functions,) to 
the boson field operator <J?(x) and the composite operator <fr(x)<fr(y). This allows the 
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determination of the one-, two-, three- and four-point connected Green's functions of the 
theory. These are manifestly finite since the underlying dynamical equations are already 
renormalized. 

Finally, we present some perspectives and discuss possible improvements of the 
approximations involved in this work. 

2 Summary of Some Exact Results 



We wish to determine in this section the exact multi-time correlation functions of 
a set of operators, which may be composite, denoted in the Schrodinger picture by Q M 
(/i = 1, 2 . . .). Following ref.|J, we introduce the generating functional 

W{£} = -ilogTrA(t )D(to), (2.1) 

where D(to) is the exact density operator of the system at some initial time to (possibly 
equal to — oo) and A(t ) is given by 



A(t ) = TexpLj™dt'J2Ut)Q»(t)Y (2.2) 



In expression fl2.2|) , £^(0 are the sources and Q^(t') are the Heisenberg representations 
of the observables Q M 

Q*(t) = U+(t,t )Q„U(t',t ), (2.3) 

U being the evolution operator. The index \i may contain both discrete indices and space 
variables that the observables Q may possess. By expanding (|2.2|) in powers of the sources, 
we obtain the following formal development for the generating functional 

+^Cd*'d* w E^U0^(O<r(Qf(0Qf(O)> + ---, 

where (O) = Tr OD (to) /Tr D (to) and = — (O) for any operator O. It is now clear 
that generates the connected Green's functions of the observables Q^. Indeed, by 

multiple differentiation with respect to the sources, one has for instance 

<Qf(0> = ^ , (2-5) 



5 2 W{i} 



(T{Q*(t)Q*{t )))=G$(t,t) = -i - (A £',,. (2-6) 

ok,ii{t )ok,u{t ) 
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The next step is to define the generating functional for the "1— PI" Green's func- 
tions, that is the effective action S e ff as a (multiple) Legendre transform of Upon 
setting 

<Q,(0> = = (Q?(0> + £~= 2 i^ 1 C dh . . . dt n 

Ei/i...^.! ^i(^l) • • • €v n -i(tn-l)G$ L ... Vn _ 1 (t ,h,.. .t n -l), 

where 

g£U(*i. • • -O = < T (Q2(*i) • • • (*»))>, (2-8) 

we can define the effective action as 

SefdiQ)} = w{0 - / dt )(Q,(t )>• (2.9) 

Now since 



SS ef f{(Q)} 



-Ut'l (2-10) 



the physical theory will be defined by the vanishing of the sources or equivalently by the 
condition 

(Q/.(0>{e = o} = <oj(0>. (2.H) 

The computation of the 1— PI Green's functions proceeds by multiple differentiation 
of S e ff with respect to (Q M (t')). But this is a rather difficult task in general since it re- 
quires the inversion of the functional relation fl2.7p . However, considerable simplifications 
occur when one intends to calculate the two-point 1— PI function . Indeed, up to first 
order in the sources, the expansion (|2.7|) yields 



at) * -* /°°dt'E (VT 1 ) (W)> - <Q?(*')>) > (2- 12 ) 

where one notes that the condition (|2.11| ) is preserved. Omitting trivial constant terms, 
the expressions for and S e ff{(Q)} now read 

W{0 ~ rdtJ2Ut)(Q^t)} + ~ rdt'dt u ^(t')G$(t',t")Ut''), (2-13) 



5 e// {(Q)} ^ 5 / dt'dt"£«Q,> - <Q?»(*)rj2(f ,0<Q*> - (Q?»(0, (2.u) 

where 



rg(t',0 s -< >k | = ( [g' 2 '1 " ) (f.O- (2.15) 



One should note that, although (|2.13|) and ( 2.1 4j ) are truncated, the expression ( 2.15Q is 
exact since is defined at the "stationary point" where the sources vanish. 
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The determination of higher order equal time Green's functions proceeds simply by 
enlarging the set of operators Q M . For instance, one may take Qi = Q (where Q is some 
simple observable) and Q2, Q3,. . . as successive powers of Q. We will see an example in 
section 6. 

The expressions that we have derived for the correlation functions require the exact 
computation of and S e ff{(Q)}. But this in turn needs the evaluation of the exact 

evolution operator which is a quite formidable task in general. In the next section, we 
present a method first developed in ref. |J and devised to provide the best approximation 
to W{£} within a variational framework. 



3 Variational Formulation 



We recall in this section the major steps in a variational derivation of the dynamics 
for the generating functional W^j^}. For a more detailed analysis, we refer the reader to 
HI sect. 3. 



The authors of || introduce the functional (with its two variants) 

l{A(t); V(t)} = TrA(t )V(t ) + 

+ j£> Tr V(t) j^ 1 - i[A®, H\ + iAit) E, Um 



(3.1) 



l{A{t); V(t)} = Tr.4(oo)£>(oo)- 



(3.2) 

In ( p.l|) and (|3.2|), T>(t) and A(t) are two time-dependent trial operators, the latter being 
subject to the constraint 

A(oo) = A(oo) = 1, (3.3) 

where A(t) is given by (|2.2j ) with t replaced by t. To account for the dynamics, the 
functional X should be made stationary with respect to variations of T>(t) and A(t). This 
leads to the following equations 



TrdV(t) 



Tr 6 A(t) 



dA(t) 
dt 

dV{t) 
dt 



-z[^(t),tf] + M(t)]T^(£)Q M 



0. 



0. 



(3.4) 



(3.5) 



One should in principle also enforce the optimization of the initial state by imposing 

TrSA(t o )[V{t )-D{t o )]=0. (3.6) 
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However, since we are not directly interested in this paper in the initial state problem, 
we will suppose that T>(t) and D(t Q ) belong to the same class so that one can impose 
V(t ) = D(t ) as initial condition. For a careful study of this problem, see || sect. 6. 

When the variations 8A(t) and 5V{t) are unrestricted, one obtains from ( |3.4j) and 

(PD 

dA(t) 



dt 
dV(t) 



- i[A(t),H] + iA(t) £ ^(t)Q M = 0, (3.7) 
+ i[H, V(t)\ - i £ £„(f)Q M Z>(*) = 0. (3.8) 



dt 



One recognizes in ( |3.7p and (|3.8|) the backward Heisenberg equation for A(t) and the 



forward Liouville-von Neumann equation for T>(t) with source terms. The stationary 
value of the action-like ( |3.1| ) reduces in this case to 

J st = Tr A(t Q )V(t ) = Tr A(t )D(t ). (3.9) 

Therefore, the generating functional (|2.1| ) writes 

W{£} = -ilogZ st . (3.10) 

Hence, one has succeeded in constructing a variational principle which has the desired 
quantity as its stationary value. 

Now, one can build approximation schemes by choosing for T>(t) and A(t) some 
trial spaces. However, the important relation ( |3.10| ) will remain true only when 5T> oc T> 



is an allowed variation [p], || because the integrand in Q3.1|) vanishes at the stationary 
point. This excludes for instance normalized operators. In addition, when 5 A oc A is 
also allowed, the quantity Tr A(t)T>(t) does not depend on the time t, as can be seen by 
adding ( |3.4|) and ( p.5|) . In this case (5V oc V, 5 A oc ^4), the approximate stationary value 
writes 

l st = expiW{£} = TrA(t)V(t) t > t , (3.11) 

as in the exact situation. Nevertheless, the dynamics will be of course different since it 
is induced by ( p.4| , |3.5|) which lead in general to non-linear, coupled equations of motion 
for A{t) and T>(t) [|, ||. Moreover, differentiation of ( p.ll|) with respect to £ yields the 
interesting result 

(3.12) 



5W{Q TrV(t)A(t)Qn 



S^(t) Tr V(t)A(t) ' 
In the following, we shall select for both T>(t) and A(t) trial spaces that meet with the 
above restrictions: the boson gaussian operators. 
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4 Approximate Evolution Equations for Gaussian 

State and Observable 



Before proceeding further, let us introduce our notation. In what follows, x denotes 
a D— dimensional vector and / x = J d D x. In addition, x is the (D + 1)— vector (x, t). 
Although the formalism presented here applies to quantum field theory, we will use a 
many-body notation, hiding as much as possible space integrals and discrete sums. 

Let us first introduce the 2N— component operator a(x) in the Schrodinger picture 

a,(x) =-^(* i (x)+*TI i (x)) j — 1,2 ... N, 

= 4| (^-at(x) - in^(x)) j = N + l,N + 2...2N, ^ 

where 3>(x) is the (N— component) boson field operator and II(x) its conjugate momen- 
tum. The usual boson commutation relations write in term of the operator a as 

N(x), «j(y)] = 7y(x, y), (4.2) 
where the 2N x 2N matrix r has the block form 

i i ( * D ( x -yA 

^) = (_^ (x _ y) o J- (43) 

We now choose for both T>(t) and A(t) the class of gaussian operators in a [ |10|] : 
V{t) = Md exp (\d,Ta) exp (iarS^a), 



(4-4) 

A(t) = M a exp (A a ra) exp (-kaTS a a). 

In (|4.4| ), Md and M a are c-numbers, A^(x, i) and A a (x, t) are 2iV— component vectors and 
Sd(x, y, i) and S^x, y, t) are 2N x 2iV matrices which can be chosen so that rSd and rS a 
are symmetric. 

The quantities M q , X q and S q (with g = d or a) uniquely define the operator to 
which they refer (T> or .4.). It is however more convenient to charecterize these operators 
by Z q (t), (a)g(x, t) and TZ q (x,y,t) defined as 

Z q (t) = TrT q (t), 

(oi) 9 (x,*) = 1 ( ; t) gV ; , (4. 5 ) 

Tr (ra)i(y)ai(x)71(t) 



(^)y(x,y,*) 
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with T d = V and T a = A. We have introduced in ( |4.5| ) the shifted operators a = a — (a). 
The relations between the two sets (Af g , X q , S q ) and (Z q , (a) q , 1Z q ) are given in the 
appendix. From now on, we consider (Z d , (a)d, Tld) and (Z a , (a) a , lZ a ) as variational 
parameters. To proceed further, we introduce the mixed operators T a d = AT) and T da = 
VA which are also gaussian operators and can therefore be characterized by (Z ad , (a) a d, 
IZad) and (Z da , (a)d a , Ttdo) given by (|4.5| ) with q = ad or da. We give in the appendix the 
relations between (Z adt d a , (at)ad,da, Tlad,da) and (Z a>d , (a) ajd , TZ a>d ). Furthermore, for any 
operator O, we will denote by (0) q (with q = a, d, ad or da) its expectation value with 
respect to A, V, T ad or T da . Our previous notation anticipated this convention. Because 
of the cyclic invariance of the trace, we have Z a d = Z da . Their common value will be 
denoted by Z. According to equation ( |3.11|) , it is precisely the quantity of interest since 
we have 

l st =expiW{Z} = Z(t). (4.6) 
The reduced functional ( |3.1|J3.2| ) now takes the form 

T = Z{to) + J Z(t) \—^-\ v + i{S ad -E da )+%Y,U^U \, (4-7) 



or 



X = Z{oo) - £° Z(t) j-^-^ - KBad ~ £*>) - * E UQ/^a j , (4-8) 

where £ a d and £ da are the expectation values of H with respect to T ad and T da 

TtHT 

e q - (H) g = (4.9) 

with q = ad or da and — \v (resp. ^~^f~U) denotes a time-differentiation in which 
the parameters of T> (resp. .4.) are kept fixed. 

The equations of motion for lZ a and (a) a are obtained from the stationarity con- 
ditions of ( |4.7|) with respect to TZ d and (a)d- A straightforward calculation yields the 
following equations 



l ^d¥ = K-n'dai 1 + Ka) - (i + K a )H ad n a , 

i^j^ = ~n a V' da + (1 + K a )V ad + t^((a) ad - (a) da ). 



(4.10) 



The stationarity conditions of ( |4.8| ) with respect to 1Z d and {a) d yield the equations of 
motion for lZ a and {a) a : 

= n d Uj\ + K d ) - (1 + K d )H' da K d , 
= -n d V ad + (l + n d )V da -t^((a) ad -(a) da ). 
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In (|4.10 , 4TTl ), we have introduced the vector V a d and the matrix 7i a d given by 



(Wad)y(x,y,t) 



dS, 



ad 



d(a) ad 

dSad. 

ad 



WZ 



(x,t) = E^i/ y ^(x,y 
(x,y,t) 



dS, 



ad 



2- 



d£ad_ 



y,ty 



(4.12) 



4J 'Wad)ji(y,X,t)" 

The primed quantities are obtained from (|4.12|) by replacing the index ad by da and the 

"energy" £ da by 

(4-13) 

This amounts to replace the Hamiltonian H by the source-dependent one H = H - 



Eu^Q,u m t ne (ia— indexed quantities as can be seen from (|4.71 , |4.8D . Furthermore, since 
Z a and Z d present no particular interest, we have prefered to omit their equations. 

The system (|4.iqj4.11| ), which to our knoweledge, has never been derived elsewhere, 
requires several comments. First of all, the coupling between T> and A occurs via the 
ad and da— indexed variables. This coupling renders the resolution of the problem quite 
complicated because of the mixed boundary conditions which are given both at the initial 
time to and at the final time t± = oo. Supposing that the initial density operator is a 
gaussian operator characterized by TZq and (a)o, we have 

K d (yL, y, t ) = 72o(x, y), (a) d (x, t ) = (ck)q(x). (4.14) 

On the other hand, the condition ( |3.3| ) implies that A/" a (oo) = 1, A a (x, oo) = and 
S a (x,y, oc) = 0. Since these two last conditions cannot be expressed conveniently in 
terms of (a) a and 1Z a (see |A.6p , it is preferable to consider the equations of motion of A a 
and T a = exp (S a ) (with T a (x, y, oo) = ^(x - y)) 

dt 



V ad -T a V da -i^(a) 



(4.15) 



I da- 



An interesting feature of (|4.10|J4.1l| or |4.15|) is that the "contractions" (7Z ad , (at) ad) 
and (TZ da , (a) da ) associated with T ad and T da satisfy the (formally) simple equations 



■An, 



ml 



A(®)ad 



[JZad, 'Had] 
Vad 



an 



d(a 



da 



t 

da 



[1Z da ,H da \ 
V' 



(4.16) 



da 



It is important to note here that these equations are not sufficient because of the mixed 
boundary conditions which cannot be implemented conveniently for (a) ad , da and lZ a d,da- 
Notice further that only (at) do and TZ da depend directly on the sources. This dissymetry 
between T ad and T da is due to the way we have introduced the source term in ( |3.1| , |3.2| ). 
Finally, upon using (|3.12|) , one obtains 

5W{{} 



(Q»)d»{t). 



(4.17) 
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This result shows that the connected Green's functions are embodied in (Q^da and not 
in (Qfj)d as one would have expected from a first view. Following the exact case, one 
may define the effective action as in section 2, with (Q M (t)) (given by ( |2.7|) ) replaced by 
(Qfj,)da(t)- The procedure now is to expand the whole set of equations in powers of the 
sources in order to determine approximate expressions for (|2.5| ), (|2.6| ) and more generally 
for the successive terms of fl2.4|) . This is the purpose of the next section. 



5 Expansion in Powers of the Sources 



When there are no sources, we have at any time A(t) = 1. Therefore, TZ a °J = 1Z 



>(o) 

da 



n w = ^(o) and ( a }W = ( a }W = ( a }W = ^(o)_ The equations (|4.iqj4.ii|gig ) become 
identical and lead to the following equations 

;dK® -[^(0) )W (0)] = _ 2 

V (o) - - ^ (0) 



■ d(a)(°) 



(5.1) 



which are the Time-Dependent Hartree-Fock-Bogoliubov (TDHFB) equations]^, [11] sub- 



mitted to the initial conditions TZS > (x, y , to) = lZo(x.,y) and (^^(x, y, t ) = («)o( x , y)- 
These equations are the bosonic counterpart of the well-known TDHF equations and 
belong to the class of the time-dependent mean-field approximation. An interesting prop- 
erty, which will be used later, is that the "von- Neumann entropy" — Tr T>(t) logP(i) is 
conserved. This leads to the fact that an initial pure state ( A.9 ) remains pure during the 
evolution. 



From expressions ( |2.5| ) and ( |4.17| ), one sees that the best variational approximation 
to (Qj(t)) within our trial spaces is given by 

(Q»}d°l(t) = (Q»}d 0) (t) = (Q^°\t). (5.2) 

Hence, to evaluate the expectation value of Q M at a time t, one runs the TDHFB equations 
(|5.1|) from to to t and then computes the quantity ( |5.2| ) by means of Wick's theorem if 
necessary. 

The next order in the sources allows the determination of the two-point connected 
Green's function Gffi. Indeed, upon using ( |2.6|) and (|4.17|) , one finds 

G^(t,t") = -z 5{Q ^ da{t '' 



5Ut 



£=0 



(5.3) 



Since (Q M )<fa depends solely on (a) da and Hda, we expand these last variables up to first 
order in the sources 



(a)da = («) (0) + 6(a) da , Tl da = K<® + 5K, 



da • 



(5.4) 
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with 



61l da (x, y, t) = -i J£° dt' E M ^(x, y, t, t%(t' 



(5.5) 



The connected Green's function GjfJ now reads 



; a(a) (0) (t) <m (0) (t) 



(5.6) 



The equations of motion for 5{a)d a and 51Zd a are obtained by developing (|4.16|) around 
((a)(°\7l^). In terms of the corrections 7 M and IZ^, they take the form 



d(a) w (ty 



z^(M') =[W M (t,0,W (p) (*)]-2 



ft(°)(t) 



(5.7) 



where 



^(t,0 = vW(*W*(m') - ^trw^^ct,/) -i<Q M >W(t)«y(t-t'). 



(5.8) 



The symbol tr in (|5.8|) , not to be confused with Tr, denotes a trace on both the discrete 
indices and the space variables. One has for instance 



2N . 



(5.9) 



The equations (\5.7j ) are not sufficient since (a) da and IZda (and hence 7 M and TZ^) are 
neither known at to nor at oo. Let us instead linearize the equations ( [4.15 ) by setting 
T a W = TP +5T a = l- S« and A« = A<°> + 5X a = 5X a , where and 5X a are of first 
order in the sources. By using a parametrization similar to ( |5.5|) , namely 



(5.10) 



L« (x, t) ee 5A a (x, t)- J y SP (x, y , t) (a) W> (y , t) 
= ~« iiT dt ' Eai £m(*')M x > t',t), 



we obtain the following equations for and 



i<^(t',t) =r 



d{a)^{ty 



AS, 



.x _ n dK,^{t ,t) 



(5.11) 



where 



MO) = V (0) (t)rL M (M) + -tr [n^(t),H (0 \t)}S,(t,t) +i(Q^°\t)8(t - t). (5.12) 
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According to the general philosophy, these equations are to be solved backward in time 
since the final conditions read 



L^(x,t, oo) = 0, S^{x,y,t ,00) 



0. 



(5.13) 



Therefore, for t' < t, we have L^(t',t) = and S^(t',t) = 0. For t' > t, L^it'^t) and 
Sfj,(t' ,t) are solutions of (|5.11|) without the 8(t — t ) term in K^. The jumps at t — t — 
read 

<9(a) (0) (x,t 

d(Q.) {0 \t 



L„(x,t ,t -0) 
S»(x,y,t',t' - 0) 



— r 



(5.14) 



We now want to derive an alternative and tractable expression for Gffi in terms of 
L M and S M . One can show first that the quantity 



g^u = L^it ,t)rj u (t,t ) - -tvSp(t ,t)TZ u (t,t ) 

does not depend on the time t except for jumps at t — t' and t = t" . Indeed, upon using 
( |5.7| . |5.11|) , one finds the expression 



^ (0)ft) 0(Q,,) u '(t) 



}S(t 



where, by choosing t" > t , it is clear that just the first term survives the integration over 
t. One recognizes in this remaining term the expression (|5.6| ) for Gffj. Hence, g^ v can be 
evaluated at t = to which yields 



G$(t',t") = L,(t',t )r lu (t ,t") - \tiS^M)n u {t Q J') 

= L„(t' ,t )TKoL„(t" ,t Q ) + £tr S^t', t )KaS„(t" , t )(l + Kq) 



(5.15) 



The second line in expression (|5.15 ) follows from the linearization of the relations between 
((a)da,T^da) and ((a) a , {a) d,1Z a ,Hd) ( |A.8|) . When evaluated at time t , the result writes 



sn da (t ) =-K s d v(t )(i + n ), 

6{a) da (t ) =KoL<p(to), 



(5.16) 



or equivalent ly 

' TZ^t") =-n S lt (t",t )(l + 'R ), 
7/n(to,t") =TZoL^(t" ,t ). 

The anti-causal Green's function is obtained from (5.15) by exchanging (/i,z/) and (t ,t 



(5.17) 
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The computation of the two-point Green's function now proceeds as follows: one 
first solves the TDHFB equations ( |5.1|) from £ to t in order to determine (a)^(t) and 
TZ^(t). Then, one runs backward ( 5.11 ) from t — t' down to t — to to obtain L M (t',to) 
and Sn(t' ,to). We note in particular that one is spared solving (|5.7|) . 

As discussed in || sect. 4, this procedure reveals several interesting features that 
are missing in the naive application of Wick's theorem. Moreover, since (Q M )da(£) and 
(Q M )(°)(t) as given by (|4.17| ) and ( |5.2|) are the best second order approximations to (Q M )(t) 
and (Q^f(t)) (respectively given by (|2/7|) and Q2.11|) ), we can use (|2.13|) and ( |2.14|) to derive 
directly the expressions 

roo j rOO 

W{£}* dt'J2Ut)(Q,) i0) (t') + - dtdt'Y,Ut)G$(t,t')Ut'), (5.18) 

Jto ,, £ Jto 



fJIV 



with the relation 



/•OO 

J to 



(5.19) 



(5.20) 



6 Illustration 



Since the TDHFB equations ( |5.1| ) remain the main ingredient whatever the set 
{Q^}, they deserve a deeper analysis. Let us focus on a 0(1) quantum field theory in 
(D + 1)— dimensional Robert son- Walker spacetime described by the metric ds 2 = dt 2 — 
a 2 (i)dx 2 . The Hamiltonian density is taken to be 

H = V D (t)II 2 (x) + a D (t)V(<!>), (6.1) 



where 



V(*) = V 2 (V<&(x)) 2 + l -{ml + g R)&(x) + j $ 4 (x), (6.2) 



with the scalar curvature given by (|| |i~2f ) R = 2D'd/a + D(D — l)(d/a) 2 . 
Upon defining the expectation values 



0(x,t) =(*(x)> 
vr(x,i) =(n(x)> 



G(x,y,t) =(*(x)#(y)> 

F(x,y,t) = (n(x)fl(y)) , (6.3) 

C(x,y,t) =(*(x)n(y) + n(y)*(x)> 
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which are related to (a) and 1Z by means of ( |A . 1 1| — [A . 1 2| ) , the energy density writes 
£ = l"~ D (t) (vr 2 (x,t) + F(x,x,t)) +a D (t)(V)(0,G). 



(6.4) 



For a gaussian density operator, (V)(<f>, G) can be computed by means of Wick's theorem 
which yields 



(V) =\{- a- 2 A x + m 2 + g R + 2A G(x, x, t) + ^0 2 (x, t)) 2 (x, t) 
+^ ( - a" 2 A x + m 2 , + + A G(x, x, t))G(x, x, t). 



In terms of the Fourier transforms of the new variables 

=a D / 2 



7T = a-^vr + ^a D / 2 



the TDHFB equations take the form 



(6.5) 



G =a D G 

p = a - D F + (Effa D G + ^-{C + C) , (6.6) 
C =C + DHa D G 



7T, 



7T 



G —C, 

F =-( a - 2 p 2 + fi 2 )C, 
C =2F -2(a- 2 p 2 + fi 2 )G, 
where we have introduced the time-dependent self-consistent mass fJ,(t): 

DH 2 , 1 



H\t) = m 2 + 



t +(g Q --)R + 2\ a- u (<p 2 + / G(k,t)), 
4 4 A 



(6.7) 



(6.8) 



(6.9) 



H = a/a being the Hubble parameter and J k = J d D k/(27i) D . In deriving ( 16.71 - 16.81 ), 
we have assumed that the field <p is homogeneous, which is compatible with translation 
invariance. 

The divergences in the k— integral of fl6.9| ) are easily removed by using the adiabatic 



expansion method[12, 13]. In the limit D = 3, this leads to the following prescription 



ni 



R _ m 



J^ = T + 77^ I T^^)■ 



l 



9r 



D- 1 
AD 



(4tt) 



9o 



D - 1 
AD 



(6.10) 



Having now established the renormalizability of the TDHFB equations, we will 
concentrate on the case where Q M= i is the 0(1) boson field operator 3>(x) and Q M= 2 is 
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the composite operator <fr(x)<fr(y). We have therefore to consider for the set a local 
source J(x, t) as well bilocal one K(~x, y, t). 

According to ( |5.10| ), L M are 2— dimensional vectors and are 2x2 matrices, which, 
for convenience, we will decompose in the form 




i ( K + u M ) 

4 \ (« M - t> M ) + 2w M 



- (u 



2tu, 



The evolution equations for these five parameters are obtained from ( |5.11| - |5"T2"D . 

read 




-d. 



,D 



( - (-a- 2 A x + M% + 2A fl 0(|C'u A 



dun 
ft 



2ia D w,, 



2ta D (4\ R 4>l tl - 
2ia D ( - a~ 2 A y 



-a- 2 A x + M 2 + 2X R g f )w lt - i\ R Cu^ 
M 2 )u„ - 2iar D v 



(6.11) 
They 

(6.12) 
(6.13) 



with M 2 = m 2 R + g^R + 2A#(0 2 + #/) and being the finite part of J k G(k, £). Their 
boundary conditions (5.14) write for /x = 1 



' ill 




,*') 


= 


ei( 


x,y,t' 


,0 


= ^ D (x 


< «i(x 




,0 


= 


ui(x 


y,z,*' 


,0 


= 


, wi(x 


y,z,*' 


,0 


= 



(6.14) 



and for fi = 2 



/ 2 (x,y,z,t 
e 2 (x,y,z,t 
w 2 (x,y,z,r,t 
D 2 (x,y,z,r,t 
. w 2 (x,y,z,r,t 





i5 D (x - y)0(z, t') + i5 D {* - z)0(y, t') 


-2(5 D (x - r)<5 D (y - z) + 6 D (x - z)5 D (y - r)) 




(6.15) 



Hence, as mentioned earlier, the optimization of the n— point connected Green's 
functions (here n = 2, 3 or 4) proceeds by solving first the TDHFB equations (|6.7H6.8|) 
forward and then, ( |6.12| - |6TT3"|) backward in time. For instance, one may solve ( |6.12| - |ErT3| ) 
with /i = 1 and \x = 2. Then, one computes G^(t',t") = G$$(x, y, t' , t"), G^(t',t") = 
G$$$(x, y, z, t',t") and G$(t',t") = G$$$$(x, y, z, r, t', t") by means of (|5.15| ). 

At equilibrium, that is for a static solution of the TDHFB equations, one can notice 
that the equations (|6.12|-|6TT3"D can be solved exactly for a wide class of metrics. One 
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can therefore check in particular that the Green's functions Gffi do only depend on the 
time-difference t — t" as expected. 

We can also notice that the two sets ( |6.12| ) and ( |6.13| ) completely decouple in the 
symmetric case = 0. Moreover, the boundary conditions (|6.14| ) imply that u±, v\ and 
wi vanish for all times t. Therefore, the two-point causal function, given by (see eq. |5.15| ) 

G$*(x,y,t',t") = -/ zr /i(x,z,t',to)i 7 o(z,r)/i(r,y,t",t ) 

-Jz,r e i( x ) z ^'^o)G'o(z,r)ei(r,y,t",to) ^ ^ 

+ 7 /z,r M x > z,t' ,t )C (z, r)ei(r, y, t" , t Q ) 
+2 Jz, r e i( x > z > t',t )C (z, T)h(r, y, t" , t ), 

where F , G and C refer to the initial density matrix 7Z , is obtained by solving just 
one equation, namely 

+ DH± + (-a- 2 A x + M 2 ))h(t, t) = 0. (6.17) 
at at 

Moreover, for fi = 2, one has only to solve ( |6.13| ) since I2 and e<i identically vanish. 



7 Conclusions and Perspectives 

We have presented in this work a variational approach suited to the optimization of 
multi-time Green's functions of a set of observables. By choosing gaussian boson operators 
for the variational objects, we obtained two sets of coupled non-linear evolution equations 
with two-time boundary conditions. The resolution of this type of equations is quite 
complicated. Fortunately, in order to compute the one and two-time Green's functions, 
one must expand these equations up to second order in a set of sources, and this in turn 
simplifies drastically the problem, since it reduces to two initial value problems. The first 
one, forward in time, turns out to be the time-dependent mean-field equations and the 
second, backward in time, is linear in the variational parameters which precisely determine 
the desired Green's functions. 

We have illustrated this formalism on a 0(1) quantum field theory in a RW metric. 
For the set of observables, we have chosen the boson field operator *fr(x) and the composite 
operator <fr(x)<fr(y). This allowed us to optimize the one-, two-, three- and four-point 
connected Green's functions of the observable 

The extension of the present formalism to fermion fields is quite simple since it only 
depends on the way we parametrize the contraction matrix 1Z (remember that in this 
case (a) = 0.) However, the generalization to gauge field theories still remains an open 
question. In particular, it is not clear whether one should introduce the gauge conditions 
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as constraints in the variational principle (which seems very natural) or one should better 
enforce these conditions in the choice of the trial spaces. The latter procedure has in fact 
been attempted [0] but has several drawbacks. 

So far we have only discussed gaussian trial operators. The question that naturally 
arises is how one can go beyond the gaussian approximation while keeping the tractability 
of the formalism and the consistency of the approximations. Some promising methods do 
indeed exist in the literature. Among them, we can quote in particular the " post-gaussian" 
approximation |4], [| which consists of an expansion in cumulants around a gaussian density 
operator, or the so-called non-gaussian calculation |~5|. More recently, a new method, 
based on the background field method, has been developed[16|]. 

We are grateful to M. Veneroni, H. Flocard and D. Vautherin for fruitful discussions. 
We are particularly indebted to C. Martin for her valuable comments and a careful reading 
of the manuscript. We would like to thank the members of the Division de Physique 
Theorique, Institut de Physique Nucleaire (IPN), Orsay- France, where part of this work 
has been done. 

In the final elaboration of this work, we received a paper [17) dealing with the same 
approach. We thank C. Martin for sending us the manuscript. 
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Appendix 



For the sake of completeness, we recall in this appendix some useful properties of 
bosonic gaussian operators of the form: 



T = A/"exp (Xra) exp (—arSa). 



(a.i; 



In ( |A.1| ), M is a c-number, A(x, t) a 2N— component vector and S"(x, y, t) is a 2N x 2N 
symplectic matrix. The (2N x 2N) matrix r is defined as 

and ct(x) is the 2N— component boson operator in the Schrddinger picture 



(A.2) 



Ct^Xj 



J_ I * i (x)+in i (x) j = l,2...N 
V2[ *j_jv(x) - zll^x) j = JV + l,...2JV 



obeying the usual commutation relations 

[«i(x), aj(y)] = 7ij(x,y). 



(A.3) 



(A.4) 



We are adopting in (|A. 1|) compact notations for discrete sums and space integrals. For 
instance 



Ara = £i ,=i / A 4 (x, t)ry (x, yK(y) 



i,j,k=l 



(A.5) 



As mentioned in the text, it is more convenient to work with the set of three param- 
eters consisting of the "partition function" Z, the vector (a) and the contraction matrix 
ft. These are given by 



Z = Tr T = A/" exp (^ArftAjVdet aft 



(a) 
ft 



Tr Ta/Z 



1 



1 — exp (-S) 



A 



A 



Tr (Tract) / Z 



exp (S 1 ) 



1 — T 

- r 

~ 1 — T 



(A.6) 



where T 



a 



a — (a) and 
o-(x,y) : 



^(x-y) 



5 D ( X -y) 



(A.7) 



The product of two gaussian operators T = T\T 2 of the form (|A.l ) with parame- 
ters Zi, (a)i, 1Z\ and Z 2 , (a)2, ft-2, is also a gaussian operator of the form ( |A.1| ) with 
parameters i?, (a) and ft. These can be written 

Z= Z 1 Z 2 exj>{^{(a) 1 - (a) 2 )rV((a) 1 - (a) 2 )}VdetaV 
(a)= KiV(a) 2 + (ft2 + l)V(a)x , (A.8) 

ft= ftiPft 2 
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with V = (1 + Tl\ + IZ2) 1 . As an immediate consequence, a projector on coherent states 
(that is a pure state density matrix): T 2 = ZT satisfies the property 

n 2 + TZ = 0. (A.9) 

The expressions of (a) and 1Z become more intuitive in terms of the usual averages 
defined as (i, j = 1, . . . N) 



&(x,t) =<*i(x)) 

7Ti(x,t) =(n i (x)) 









)*i(y)> 






= (n^x 


)n,(y)) 




&y,t) 




in^ + n^y^x)) 



(A.10) 



with Q = Q — (Q). Indeed, with the help of the canonical transformations ( |A.3| ), one gets 
easily 

(Q '^7! (A ' U) 

and 
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